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Abstract. Starting from the generalized t — J — G model Hamiltonian, we analyze the spin response in the 
superconducting cuprates taking into account both local and itinerant spin components which are coupled 
to each other self-consistently. We demonstrate that derived expression reproduces the basic observations 
of neutron scattering data in Y Ba2CuzO(,+y compounds near the optimal doping level. 
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1 Introduction 

Theoretical studies of layered cuprates can be broadly sep- 
arated into two parts depending on whether the point of 
consideration is Mott insulator or a metal. Theories of 
Mott insulatorsfl] departs from a Mott insulator/Heisen- 
berg antiferromagnet (AFM) at half-filling, and addresses 
the issue how superconductivity (SC) and metallicity arise 
upon doping. Another class of theories [2 explores the 
idea that the system's behaviour is primarily governed by 
interactions at energies smaller than the fermionic band- 
width, while contributions from higher energies account 
only for the renormalizations of the input parameters. 
For the cuprates, the point of departure for such theo- 
ries is a Fermi liquid (FL) at large doping, and the issue 
these theories address is how non- Fermi liquid physics and 
unconventional pairing arise upon reducing the doping. 
Correspondingly, two approaches are also used for the de- 
scription of a dynamic spin susceptibility. In the itinerant 
case usually one employs the random phase approximation 
(RPA) 3,4,5,6,7][5] and in a proximity to the antifcrro- 
magnetic state, the spin response above Tc is governed by 
the continuum of the antiferromagnetic spin fiuctuations 
(paramagnons). In the superconducting dx2_y2-wsiVQ state 
there is a feedback effect of superconductivity on the spin 
response which yields a formation of the spin resonance 
below Tc. While the behaviour of the spin response in the 
superconducting state of optimally and overdoped cuprate 
superconductors can be qualitatively and even quantita- 
tively understood within this approach, the normal state 
data are not entirely captured by the RPA where the ex- 
citations are completely damped and structureless at high 
energies. 
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Within the localised type of approaches the situation 
is opposite. In this case one starts from the undoped sit- 
uation of a two-dimensional antiferromagnet and studies 
how the spin excitations evolve upon introducing the finite 
amount of carriers 9, 10, 11 , 12, 13, 14, 15, 16^. Here the spin 
response remains dual in nature as it assumes a mixture 
of the local spins described by the superexchange interac- 
tion J and the itinerant carriers with tight-binding energy 
dispersion. This scenario seems to be more efficient in de- 
scribing the normal state spin dynamics, but so far its 
application to the superconducting state was rather lim- 
ited. So called downward dispersion of neutron scattering 
intensity is not reproduced in this approach. 

Note, the hour-glass-shape dispersion observed in neu- 
tron scattering below Tc naturally calls for the explana- 
tion of the spin response in terms of dual character of 
the excitations. While the upward dispersion resembles 
the collective spin wave-like branch as in quasi-two di- 
mensional antiferromagnet with short range spin fluctu- 
ation, the downward dispersion in the superconducting 
state can be nicely attributed to the feedback effects of 
the d-wave order parameter on the itinerant component. 
However, usually an interaction between both would in- 
troduce the repulsion between both branches and it is not 
'a-priori' clear how the total response would look in this 
case. In this paper we discuss the possible way how to 
describe both components (local and itinerant) on equal 
footing within one analytical scheme based on the Green's 
function method. 



2 Basic equation for quasiparticle operators. 

The starting point for our analysis is the usual t — J type 
Hamiltonian with additional density-density interaction 
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term 



I J I Si Sj 



(1) 



Here, if^P'^-'^^if-j^'P^) are the creation (annihilation) oper- 
ators for the composite quasiparticles within the conduc- 
tion band of hole-doped cuprates written in terms of pro- 
jective operators to obey the no-double occupancy con- 
straint. The second and third terms describe the superex- 
change interaction between the spins and the screened 
Coulomb interaction of the doped carriers, respectively. 
(5 is a number of holes per one unit cell. Note that the 
spin operator commutes with the density-density inter- 
action and this it will not appear explicitly in the final 
expression for the spin susceptibility. 

Let us begin with the equation of motion for Fourier- 
transform of the quasiparticle operator 



(2) 



The linearization of the commutator on the rhs of Eq. 
^ can be performed via projection method on the space 
of creation and annihilation operators, in a similar way 
proposed previously fl7','18''19I. However, in contrast to 
the previous approaches [1 2, 16| we keep also the molec- 
ular field terms, which are proportional to the Fourier- 
transform of the transverse spin and the density of holes. 
In particular, we approximate the commutator in ^ by: 



Note, that the energy dispersion, Eq.Q, can be writ- 
ten in a conventional form of the tight-binding dispersion. 



£k = 2t|. J J (cos k^a + cos kyO) + 'i^e'f'f cos fc^a cos kyU. (5) 
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but with effective hopping integrals 
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where for shortness we introduce P = and the spin- 
spin correlation functions Kn = 4 {S^S^) are calculated 
via the spin susceptibility expression in a self-consistent 
way for a given doping and temperature. 

In the following section we will further proceed with 
the derivation of the dynamical spin susceptibility based 
on the basic equations of motion for the quasiparticle 
operators. This approach differs from the memory func- 
tion method (MFF) [l2l[T3l[T6] where a linearization of 
the equations for the Fourier-transform of the spin oper- 
ators is not related to the basic equation for quasiparti- 
cale operator ([3]). In this regard our approach is closer to 
the conventional random phase approximation. However, 
it employs projecting operators which fulfil the anticom- 
mutator relation 



1 



(8) 



dt 



N 



q^k-q ^q 



1 -(T^ , 

+ 2iV^** 



k-qVf,_gdq 



2N Z^'^l^k-q^q 



+ 



-T 



Gn 



^k~q "q' 



(3) 



where t'^. = tij{l — Fj) exp {ikRij), Jq — 2Ji{cosqxa + 
cos qyu) and Gq = J2j ^ij exp(iq'i?y ) are the Fourier-trans- 
forms of the hoping integrals, superexchange, and Coulomb 
interactions, respectively. In the following we also set the 
lattice constant a to unity. is a number of Cu-sites 
in the copper-oxygen plane. Collecting terms which are 
proportional the hopping amplitude, tji , the following ex- 
pression for the energy of quasiparticles is obtained 



Sk 



E' 



l + S 



l + S 



(4) 



The physical meaning of square brackets can be under- 
stood as follows. The antiferromagnetic spin correlations 
suppress effective hopping integral, while ferromagnetic 
ones increase it. Thus one should expect that 1 + 2Fi > 0. 
Note that overall the effective nearest neighbour hopping 



determined as 



l+S 
2 
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1+.5 



(1 + 2Fi) {S^Sf] 



to Gutzwiller's projection for the hopping integral 
[20] . We show below that the exact form of Fi can be 
rigorously computed via carriers concentration S and the 
spin-spin correlations functions, see Eg. (1281) . 



is analogous 



which results in richer behavior of the spin susceptibility. 
Although this approach requires some intuitive knowledge 
on the behaviour of the physical system it has some advan- 
tages because (i) both static and dynamic susceptibilities 
can be calculated within one approximation scheme, based 
on ([3]) , (ii) the extension of the formalism to the various 
symmetry-broken ground states such as superconducting 
or spin/charge density wave ordered states is straightfor- 
ward by adding the corresponding terms in the linearized 
equation of motion, Eq.(3), and (iii) the correct asymp- 
totic Fermi-liquid type behaviour of the susceptibilities, 
when spin-spin correlation functions become small and the 
system approaches the conventional Fermi-liquid regime is 
restored. 

In particular, Eq.(3) contains the superconducting gap, 
Z\k which mean-field expression can be readily found 



— T 

PN ^ 



Jk-k' i^ll^'^ttf) 



PN 



E 



-k' 



-tl. 



k' 



(9) 



where tk' = J2j *u ex.p{ikRij), t^, = J^j UjFj eyip{ikRij) 
and = Gq — Jqj'l- One sees that both superexchange 
and density-density interactions are involved. The analy- 
sis of this equation for Jfe_fe' > Gk-k' , which we assume 
hereafter, reveals dj.2_j,2— wave symmetry to be the most 
stable solution with Ak = A(T){coskxa — coskyo) /2. 
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3 Dynamic spin susceptibility. Combining Eqs.(ITll)-([T7]) we find 

A hierarcliy of constructed equations of motion for deter- d'{oj, q) = ~ \ <?) + q){{St\SZ„)) \ K{u], g) 



mining the spin response contains five Green's functions. 
The initial one in the absence of the long range order 
{{S,)^ 0) has the form ^'^^ 



(20) 



k' 

+ E J^.le-''"'' {{S^S! - S!St\SZ,)), (10) 



where tk is the usual tight-binding Fourier-transform hop- 



(21) 

where for T > 



ping integral on the square lattice including the nearest, i i 

next-nearest, and next-next-nearest neighbour hoppings C(^.g) = / ,Ck.q ^ Tr / . ^ | ' 

(a bare dispersion). The second term 

Gioc{u;,q) = E Aie-^^'HiS^S^ ~ S[-S+\SZg)) (11) v'{^,q) = E'?^ = ^J^xi^^l) 

k 

refers to the contribution of the localized spins. Its form ^k+gi^k+g — tki^k ^^^-j 

is determined by the same procedure as described in Refs. N ^ lo ^ £k — Sk+q 

13110,11,12,15^ 



<^GiU^,q) = ~-Jiifi(2 - 7,) + nl{{S+\SZ^)) (12) 

TT 



V {^,1) ='^V'k.q = V'i^,<l) 



and Hq determines the frequency of collective local spin + — ^(^k+q — *fe) — ^ /25'j 

' ■ / 1-1 \ N UJ + £k-£k + q 



fluctuations (magnon-like) 

= 2J?q(2 - 7,)(zi,p -K,{2 + 7,)), (13) 



We remind that the difference between tk and ij, is ex- 
pressed below Eq.(3) and Uk = Pfk where fk is a Fermi 



where Z\,p is a dimensionless parameter of the so-called faction. Takmg mto account that Efc' i^k'+q - Efe') = 

spin-gap, [T21IT3] , 7, = cosg^a -h cosg^a, and a is a de- it is easy to prove that D' (w, q) = D (w, q). This means 

coupling parameter, controlled by the sum rule {S+Sr) = that Eqs.iEOl) and (EH) are consistent with each other if 

1 /, i-N . the following relation holds 

^(1 — 0), which IS about 1.4. ° 

The first term on the right hand side of Eq. ((10]) 1 P{t'k+q - t'k) 



(26) 



2 . ... 

k' This property can be used to express the projection pa- 

is determined by the dynamics of itinerant spins. To com- rameter Fi via the spin-spin correlation function Ki. For 
pute this function we use the following exact relations [HI this we add and subtract Sk - Sk+g m the numerator of 
[52] Eq.(26) which yields the identity 

j:^fMt:-0, (14) ^ P{t'k^q-ti)~Ek^q+Ek p 

* ^ LJ + Ek- £k + q 

^^t.P.^P.4^5+ (15) ^ , 

fe Note that the remaining sum vanishes at any frequency 

Using Eq. m we construct further equations for the ^^^'^+1 " " ^^+« + = 0. Equating factors in front 
Green's functions '^^ independent trigonometric cosine functions results 

in the relation between Fi and corresponding spin-spin 
m'-'^lfq\SZg)) = l-Xk.q+LnUiSt\SZ,)) + ^-fD' {u,q) , correlation functions 

= -^;,,,^_^,;'_,((s+|5Z,»-%iD" (., q). (' + + '''' 

^lY) which holds for any value of i. Recalling that 1 -|- 2Fi > 
Observe that here the new Green's functions appear 0, we find that (1 -f- > 2\Ki\. This condition implies 

that the concentration of carriers should be larger than 

D' (a;, g) = - E ('^fc'+9 ~ efc')((*'fc''*'^*'fe'+'!,l'S'-q)). (18) some critical value Sq to have the metallic behaviour of 

k' the system, 

n" „\-Sr(^ P M/m^'P'^mP'i'i-lQ- \\ ^^Q^ Having D' {uj,q), D" {uj,q) and using Eqs.dH]), ^ it 

^ V " " '^^ " ' ^ ^ is straightforward to find Gu{uJ,q) and ((5+|5Z,)). The 
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obtained system of coupled equations allow us to get a 
closed form for the Green's function {{S'^\SZq)) i.e., for 
the spin susceptibility 

+,- _ ^x('-^:g)Ct(t^,g) - [1^x^(1^. g) + 2/1-^1(2 -7g)]c(a;,g) 



uj[l/2 + Tj{uj,q)](t{io,q) + - uj^ - a;r?t(tj, g)]C(a;, g) ' 



Here, 
and 



-JgXii^,q) 



J2 - -P/2) - - P/2) 

k 



IjJ + Ek — £k+q 



(29) 



(30) 



Xt{oJ,q) 
Tlt{oj,q) 



(31) 



- tk)rik,q. 



Note that the analytic continuation lo = lo + iO^ used 
in the derivation is replaced in the numerical calculations 
by w = ui + ir where -T is a small numerical factor. It 
is frequency independent and resembles the effect of the 
non-magnetic impurity scattering. [51I21[Z1[S] 

Let us turn now to the situation of the superconduct- 
ing state, T < Tc. Performing the Bogolyubov transfor- 
mations for the quasiparticle states and obtaining the new 
Eigenenergies Ek, we then employ the new projecting op- 
erators, that takes into account the {it, v} Bogolyubov 
coefficients. This procedure is standard in the theory of 
superconductivity and thus we skip the details. Note that 
overall expression for the spin response retains its form 
except that the entering functions such as x(a;,q) modify: 



TV 



fk 



fk + q — fk 

uj + Ek — Ek^ 

fk + Q 



LJ ~ Ek + Ek + q 
fk + q + /fe — 1 



+ 



(32) 



Ek — Ek + q 

-A + ) 1 ^ fk + q + fk 

LU + Ek + Ek + q ' 

which resembles BCS-like expression for the spin suscep- 
tibility for non-interaction electrons and the Fermi func- 
tion fk contains now the new eigenenergies. Note that the 
superconducting gap is determined by the numerical so- 
lution of Eq.(9). For the sake of simplicity, we further use 
the following abbreviations for the Bogolyubov coherence 
factors: 



where 



Sxx 

(+) _ 



- XkXk + q ~\~ ■^k^k + q: 
XkVk+q — ZkZk+q, 



ykVk + q + ZkZk + q, 
TJkXk + q — ZkZk + q, 



(33) 



Zk 
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Ak 
2Ek 
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Vk 



Ek 

Ek = V{£k-fir + \Ak 



Ek 



(34) 



Then the function i]{uj,q) gets the form 
1 P 

Sr c t'k+,{fk+q- 1/2) -t'kifk- 1/2) , 

^ UJ + Ek- Ek + q 

^V- . t'k^,{l/2^fk+q)-t'k(l/2~fk) , 
+ ™ co-Ek + Ek+q + 

, ^(-) 4+,(/fe+.-l/2)-4(l/2-/0 
^ V uj-Ek- Ek+q 

.sr^ c.(+) 4+g(f/2 - fk+q) - t'kifk - 1/2) 1 

^ UJ + lijk + Ek + q J 

where P is again the thermodynamic average of anticom- 
niutator 

in the normal state. As mentioned above the anticommu- 
tation relations for projecting operators are not simple 
fermionic ones. As a result the last terms which appear 
in Eq. ([55]) correspond to an effective molecular field of 
kinematic origin which is a result of the projection origin 
of the operators (strong correlation effects) . 
The function C(w,q) is written as follows: 



N ^ 



N ^ u + Ek-Ek+q 



N ^ UJ-Ek + Ek + q 



1 

^ N ^ uj-Ek-Ek'^ 



k — J^k + q 



+ ^E 



s'. 



(+) 



U + Ek+Ek + q 



(36) 



At the same time the terms in Eq. ([25] which describe the 
coupling between the magnetizations of the itinerant and 
localized spins remains the same as in the normal state, 
see dSU. Observe also that Eqs. dSS]), dSp refer to 

itinerant spin-component of the spin susceptibility. If we 
assume for the moment that no hopping terms ( i.e.ii = 
t2 = 0) of conduction electrons (holes) are present, Eq. 
((^ reduces to 



Xloc 



2JiKi{2^-yq) 



(37) 



This expression is identical to those found previously and 
is widely used to describe the hghtly doped cuprates [TOl 
[TT1[T^[T51[T3] . It is remarkable that magnetism of localized 
spins in Eq. (|29p is strongly suppressed (or in other words 
"frozen out") due to the superconducting gap, which is 
naturally incorporated in expression for ({uj,q). In the 
opposite limit, when the spin-spin correlation function is 
small and conduction electrons bandwidth is large enough 
(C(a;, q) is small), Eq. (1251) becomes similar (but not identi- 
cal) to the one obtained previously in the so-called gener- 
alized random phase approximation (GRPA) scheme for 
Hubbard model in the normal state [23] • In contrast to 
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the previously obtained expression, see Ref. [33], Eq.(29) 
obeys the electron-hole symmetry which is especially im- 
portant for the superconducting state. 

In addition note that Eg. ([3^ contains contributions 
from both the itinerant and the local components of the 
spin susceptibility, which are mutually coupled in a non- 
trivial way. In fact, these components cannot be easily sep- 
arated and should be treated as collective spin response 
of the entire system. The energy position of the spin exci- 
tations is obtained by analysing the pole structure of the 
denominator of Ea. (j29p . In quasi- localised spins regime 
this equation corresponds to a short-range magnon-type 
oscillations with an upward dispersion near the antiferro- 
magnetic wave vector Q = (tt, tt). In the case of the itin- 
erant spins the expression (j29l) yields a Stoner continuum 
with the overdamped paramagnon excitation. 



4 Numerical results and comparison to 
experimental data 

In the following we present the numerical results using 
Eq. ([35)1 for the normal and superconducting state. As 
we pointed out above the superconducting gap equations 
yields 3,^2 _y2-'wsNe symmetry of the superconducting gap 
for Jk-k' > Gk-k' ■ Therefore we approximate the super- 
conducting gap in the form Ak = (cos k^a — cos kya) /2. 
The magnitude of Aq w 30mey is obtained from the 
temperature dependencies of the nuclear relaxation rate 
[23] and superfluid density [53] for YBa2Cu30r. The en- 
ergy dispersion is given by (jSJ and we employ the follow- 
ing minimal set of effective hoping parameters (in meV): 
tg^j = 250, tgj-'y. = —50 which reproduce the observed 
Fermi surface for optimally-doped cuprates near the opti- 
mal doping level. 

The calculated imaginary part of susceptibility in the 
normal phase for T — lOOK is shown in Fig. \T\ The 
chosen parameters arc: Ji = lOOmeV, A^p = 0.18, and 
r = 6meV. Values Fi = 0.15 and F2 = -0.02 was es- 
timated using Eq (gS]) at Ki = -0.2, K2 = 0.04 and 
6 = 0.3. Observe that the imaginary part of total suscep- 
tibility has a visible upward dispersion, which can be at- 
tributed to a magnon-like mode damped at high frequen- 
cies due to coupling to the itinerant carriers. Note that 
the visible dispersion of the spin excitations centered at 
the antiferromagnetic wave vector Q, is in contrast to that 
found within a simple RPA expression [2 6^. There the spin 
excitations are almost structureless in the normal state 
and refer to the overdamped Stoner continuum. Another 
important feature is that the spin excitations remain com- 
mensurate despite the significant hole doping as shown in 
Fig[TJa). This is due to the fact that the spin excitations 
has a character of the almost localised magnetic modes 
which are less sensitive to the degree of the Fermi surface 
nesting. One has, however, to keep in mind that in our cal- 
culations we take 7^ = 6meV as a constant. In principle 
in a real physical system this value may become momen- 
tum and frequency dependent. In particular, assuming the 
coupling of the quasiparticles to the spin excitations, the 



quaiparticle lifetime will become anisotropic as a function 
of the Fermi surface angle. Quasiparticles connected by 
the antiferromagnetic momentum will scatter stronger as 
compared to those located near the diagonal of the BZ. 
As a result, the imaginary part of the self-energy (which 
affects r) will be larger around the wave vector Q and is 
smaller away from it. To demonstrate this effect we ap- 
proximate P by a Lorenzian in the form 



r{q) = A + 1^2 



i\q.\-^r + i\Qy\~^r + iO' 



(38) 



where ^ is a magnetic correlation length. As one sees 
the effect of q-dependence of the damping F yields a 
dip in imaginary part of susceptibility at (tt, tt) as shown 
in FiglTJb) and to the incommensurate magnetic fluctua- 
tions. 

In Fig. [5] we show the imaginary part of the spin sus- 
ceptibility in the superconducting state with d^^ -y2-waye 
symmetry. Observer the strong renormalization of the en- 
tire spectrum of the spin excitations. As a consequence 
of the d-wave symmetry of the order parameter there is 
a resonance mode forming at energies well below |Z\k -I- 
^k+ql ~ 60meV. Away from the antiferromagnetic mo- 
mentum, the excitations shows the X— shape structure 
which appears only in the superconducting state. Due to 
the delicate balance between the spin-gap parameter Ag 
and the twice superconducting gap, 2Aq, the frequency 
position of the resonance peak in superconducting state is 
almost unchanged with regard to the normal state value. 
However, its position is still below the onset of particle- 
hole continuum which starts approximately at 2Ao. Tak- 
ing the constant energy cuts of the excitations shown in 
Figs. [T] and one finds that the results are very simi- 
lar to those observed in inelastic neutron scattering in 
YBa2Cu30Q,g2 and YBa2Cu30(i,g7 (see Fig. 8) in P. Bourges 
Rev. _27 . Note that the calculated upward and down- 
ward dispersions shown in Fig. [5] look similar to those 
found experimentally in Y Ba2Cu^O(, 95 (see Fig. 3(c) in 
M). (Fig. 5(a)-(b) in Ref. [H]), FBaaCT/aOe.ss (Fig. 
4(a) in Ref. [25]), YBa2Cu30e.6 (Fig.5 in Ref. [HU]), and 
YBa2Cu30Q,5 (Fig. 5 in Ref.f31j'). The coexistence of both 
downward and upward branches can be clearly seen in 
Fig. [2] While the downward dispersion refers to the itin- 
erant component of the spin response and the evolution of 
d-wave gap on the Fermi surface, the upward dispersion 
originates mostly from magnon-like short-range antiferro- 
magnetic fiuctuations. Important to note here is that these 
dual features are obtained using the general expression for 
the spin susceptibility which accounts for both local and 
itinerant spin excitations. 

In Fig 121 we show the intensity plot of the imaginary 
part of the susceptibility /mx+'~(tj, q) as a function of q 
around the wave vector Q = (tt, tt) in the superconduct- 
ing state for the energies above and below the resonance 
energy positions, w^es- Below w^es in agreement with ex- 
perimental observation and RPA calculations IIT we 
have found quadratic in q downward dispersion with dom- 
inant peaks along the bond directions. The latter is due to 
fact that there is more phase space available for (tt, ti + S) 
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Fig. 1. Calculated imaginary part of the spin susceptibility for the normal state (T=100K) near the wave vector Q — {n,n) as 
a function of a; (in eV) and Qx [q-y = n] for F — 6mcV (a) and that using Ea. H38ll with A — 6meV, J2 ~ QmeV , ^ = O.OStt. (b) 




Fig. 2. Calculated imaginary part of spin susceptibility for the superconducting state (T=10K) near the wave vector Q — {n,n) 
as a function of to (in eV) and for qy = n. 



directions as they are closer to the resonance at (tt, tt) 
than those at the diagonal wave vector at this energy. At 
the same time, the upward branch above w^es has almost 
circular symmetry around Qaf which is the same as in 
the normal state and agrees with experimental observa- 
tion [26j . In order to see the evolution of the incommen- 
surability with the superconducting gap we plot Aq vs 
26 /tt in Fig. m Observe that for increasing Aq the incom- 
mensurability parameter, S, computed for uj = AOmeV in- 
creases. This is an indication that the incommensurability 
for UJ < LOres rcfcrs to the spin exciton dispersion which is 



connected to the magnitude of the superconducting gap. 
Namely, the increase of Aq shifts the position of the res- 
onance peak to higher energies. For a given w < ujres this 
increase implies that the itinerant spin excitations arises 
from the scattering of the quasiparticle states located at 
the Fermi surface points which are closer and closer to the 
diagonals of the BZ, i.e. smaller (larger 2^). Note that 
this observation is in direct correspondence with the result 
found previously in experimental data (see Fig. 25 in Ref. 
|32|). At certain point, when the resonance energy is large 
enough, the only electron-hole scattering available orig- 
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Fig. 3. Intensity plot of the imaginary part of tiie susceptibility 
Imx'^'~ iijJ,q) around wave vector Q = {tt,tt) for huj — 40meV 
(a) and hw = QOmeV (b),r = 6meV. 



inates from the nodal points where the superconducting 
gap is zero. At this moment, the incommensurability starts 
to be independent on the Aq, the effect clearly visible in 
FigSl This observation confirms the itinerant origin of 
the lower branch of the dispersion of the spin excitations. 
Note that the incommensurability of the spin excitations 
for CO > LOres IS uot Sensitive to the magnitude of the su- 
perconducting gap as it arises from the localised magnetic 
excitations. In this case it remains a constant (not shown) 
independent on the size of the superconducting gap. 

In Fig|5] we show the relation between the resonance 
energy and superconductivity gap Aq. Despite the fact 
that this observation is a direct consequence of the reso- 
nance being the spin exciton we show it here as it agrees 
fairly well with experimental data (see, for example. Fig. 
26 in Ref. [21]). 

Finally, in Fig|6]we show the calculated temperature 
dependence of the resonance peak using the temperature 
dependence of the gap 




25/7C 

Fig. 4. The incommensurability parameter, 25/n plotted as 
a function of the superconducting gap amplitude Ao- The red 
points refer to the calculated values using Eq. p9|l and the 
solid line is a guide to the eye. 
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Fig. 5. Calculated relationship between the resonance energy 
(in meV) and superconductivity gap zio- 



A{T) = Aot8inh(l.76^/TJT~^^ which was obtained 

Dreviouslv [24ll25| . The resonance peak intensity follows 
the temperature dependence of the superconducting gap. 
Above Tc the resonance mode continuously evolves into 
paramagnon mode. We stress that the true resonance ex- 
ists only below Tc and the excitations above Tc refers to 
the spectrum of the normal state. This result is in agree- 
ment with the experimental observation found in Y Ba2Cu^OQ,g'j 
(see Fig.lO in Ref. [27j and also Fig.4(a) in Ref. [33]). 

At the end of this section we shortly summarize the 
differences between our results and those obtained earlier 
either in a conventional RPA scheme or within a mem- 
ory functions formalism (MFF). As it was pointed out in 
previous RPA studies, see Ref. [26], the conventional RPA 
scheme explains quite well the downward part of the dis- 
persion of the spin excitation which originates in this case 
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Fig. 6. Temperature dependence of the resonance peak inten- 
sity. 



due to d-wave symmetry of the superconducting gap, but 
does not reproduce the upward one which should refer to 
the paramagnon dispersion. The origin of this is that the 
lifetime of the quasiparticles is not included within sim- 
ple RPA. As a result, the paramagnons at higher energies 
are strongly damped due to a standard Landau damp- 
ing mechanism. In our approach the Landau damping is 
suppressed by the kinematic " molecular" field which origi- 
nates from local no-double occupancy constraint. The con- 
straint reduces the phase space available for the particle- 
hole continuum to form and suppresses then the Landau 
damping. As a result the upper paramagnon part of the 
dispersion is more visible. 

At the same time the chains of equations of motion 
with a decoupling procedure employed by us and mem- 
ory function formalism used by other groups are closely 
related. In the memory function formalism one decouples 
the exchange part of the Hamiltonian in fashion similar to 
ours. The resulting dispersion of the localized excitations 
is then again referred to the upward branch of the spin 
excitation seen in INS. To account for the superconduc- 
tivity and presence of itinerant electrons, which give rise to 
the downward dispersion, one introduces the decoupling of 
the spin-excitation damping. This describes the decay of 
the spin excitation into an electron-hole pair. Here, super- 
conductivity is introduced phenomenologically in the en- 
ergy spectrum of the quasiparticles. One finds that the re- 
sults for the spin excitations within memory function and 
equations of motion formalism are quite similar, see Refs. 
P^IT^ITC] . However, we believe that our approach dealing 
with the chains of equations of motion has certain advan- 
tages. In particular, it allows to reproduce wider range of 
experimental data on INS for T <Tc (also including NMR 
and ARPES data) and, most importantly, we are able to 
describe both upward and downward dispersions (see Fig. 
3) within a single analytical expression. Here, the origin 
of various terms have a transparent physics meaning and 



the parameters of the quasiparticle dispersion and the spin 
excitation spectrum are computed self-consistently. 



5 Concluding remarks 

To conclude, in this paper we present the analysis of the 
spin response in the superconducting cuprates taking into 
account both the local and the itinerant component of 
the magnetic excitations coupled in a self-consistent man- 
ner. The numerical results show that the obtained ex- 
pression for the spin susceptibility, Eq. (f^ . reproduces 
well the characteristic features of the experimental data 
in Y Ba2Cu^0Q+y compounds near the optimal doping in- 
cluding the dispersion of the spin excitations in the normal 
and superconducting state as well its frequency and tem- 
perature dependence. While the structure of the spin exci- 
tations in the normal state can be attributed to the over- 
damped localised magnetic modes, the strong feedback of 
the d-wave superconductivity on the itinerant electrons 
reveals the formation of the spin exciton with the charac- 
teristic downward dispersion of the spin excitations in the 
superconducting state. In addition, the high energy spin 
excitations still originates from the localised excitations. 
Remarkable that both modes merge at the uJres which is 
a result of the single pole structure at Q in Eq. ([^ . Our 
analytical results suggests the dual character of the spin 
response even in the optimally doped cuprates. 
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marks and discussions. This work was supported in part 
by the Russian Foundation for Basic Research, Grant 09- 
02-00777-a and Swiss National Science Foundation, Grant 
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